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Abstract 

We study the dependence of the phase behavior of ternary amphiphilic systems 
on composition and temperature. Our analysis is based on a curvature elastic 
model of the amphiphile film with sufficiently large spontaneous curvature cq and 
sufficiently negative saddle-splay modulus R that the stable phases are the lamellar 
phase and a droplet microemulsion. In addition to the curvature energy, we consider 
the contributions to the free energy of the long-ranged van der Waals interaction 
and of the undulation modes. We find that for bending rigidities of order /c#T, 
the lamellar phase extends further and further into the water apex of the phase 
diagram as the phase inversion temperature is approached, in good agreement with 
experimental results. 
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1 Introduction 



The theory of ternary amphiphilic systems has developed rapidly during the last few 
years Jl], [2|]. For long-chain amphiphiles the molecular solubility is very small, so that 
essentially all amphiphilic molecules are located in monolayers at the interface between 
oil and water. In this long-chain limit, the description of the system by an ensemble of 
surfaces is justified and used frequently 0. The monolayer is described by an ideal surface 
in this approach, with the shape and fluctuations determined by the curvature energy |3J] 



where H and K are the mean and Gaussian curvatures, respectively, and the integral runs 
over the total interface area. The parameters are the bending rigidity k, the saddle-splay 
modulus R, and the spontaneous curvature cq. 

Most of the previous theoretical work on the phase behavior of ternary amphiphilic 
systems has focused on balanced systems with Co = 0, and on the interplay between the 
curvature energy, which favors the lamellar phase, and entropy, which favors disordered 
bicontinuous or droplet phases |L], [|. However, in addition to the curvature energy, 
long-range interactions between the molecules also give an important contribution to the 
free energy. This is well known for the unbinding of membranes ||, and for wetting 
transitions in simple fluids || . Attractive interactions have also been used to explain the 
reentrant phase separation of spherical and (very long and flexible) cylindrical micelles 
near the water apex of the composition phase diagram Jj]]. In systems containing non- 
ionic amphiphiles, or in systems with screened electrostatic interactions, the dominant 
contribution is the van der Waals interaction || |^]. 

In this paper, we want to study the dependence of phase diagrams on the spontaneous 
curvature c . We consider systems with c sufficiently large that the isotropic phase is 
a droplet microemulsion (in order to avoid complications with the bicontinuous phase). 
It has been shown experimentally that cq varies linearly with temperature T over a wide 
range of T flOfl . Using this relation for co(T), we can compare our results directly with 
the temperature dependence of experimental phase diagrams. 




(1) 



2 



It has been shown [ITJ that in the zero-temperature limit, the curvature model - 
together with the constraints of fixed volume fractions <p s of the amphiphile and <p of 
the interior phase — predicts for —0.5 < k/k < with increasing Co<j) /<j) s the progres- 
sion from a lamellar phase via a phase of infinitely long cylindrical micelles to a phase 
of spherical micelles, while for k/k < —0.5 the phase of cylindrical micelles is absent. 
This model has been extended to include both thermal fluctuations and long-range in- 



teractions in Ref. [12]. However, with the interactions described phenomenologically, and 
no explicit temperature dependence, this model yields highly schematic phase diagrams, 
which cannot be compared easily with experiment. 



2 The Model 

We consider the coexistence of a lamellar phase L a and an oil-in-water microemulsion 
phase L\. We assume that R is sufficiently negative that a hexagonal phase of cylindrical 
micelles is suppressed [IT]. The L^-phase is a stack of oil- and water layers of thickness 



d and I, respectively. All lengths are measured in units of the amphiphile length, so that 
the repeat distance in the lamellar phase is given by L = d + I + 2. The radius of the 
neutral surface of the oil droplets is denoted R, so that the radius of the oil core is R— 1/2, 
while the outer radius of the micelle is R + 1/2 (for equal lengths of heads and tails). The 
number density of droplets is p. 

We write the free-energy density, /, of the lamellar and the droplet phases as a sum 
of three contributions, 

f fcurv f therm fvdW- (2) 

The curvature-energy density, fcurv, is given by 

I 4/tCo/L for L a 

fcurv \ (3) 

[ 4vrp (2k(1 - c R) 2 + R) for L x 

Thermal fluctuations in the lamellar phase lead to a steric (repulsive) interaction, as first 
proposed by Helfrich [|13| . When the size of the amphiphile is used as an ultraviolet cutoff 
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for the wavelengths of the undulation modes, the steric-interaction density of pairs of 
monolayers, which are separated by a water film of thickness I, is found to be 

UncM = ^ [ti Ml + (^||)- 4 ] + | arctan[(g c £||) 2 ] \ , (4) 

where q c = 2n, and £n = (Sup/hsT) 1 ! 2 I is the parallel correlation length. The parameter 
p has been determined by Monte Carlo simulations to be p ~ 0.135 JT1, |T5]. Note that 
Eq. (|j) reduces to the familiar / _2 -form of the steric interaction of two monolayers 
for qr c £|| > 1. 

It has been pointed out recently that for an asymmetric lamellar phase, the steric 
interaction of two monolayers is a function of both the thickness of the oil and of the 
water films W&. This is obvious in the case of vanishing oil concentration, where the 
interacting surfaces are amphiphilic bilayers with an effective bending rigidity K e ff = 2k. 



In Ref. ||16|| , Monte Carlo simulations have been used to extract the expression 

K eff (l; d)/n = 2 - (d/iy/ [(d/iy + x v ] (5) 

with v = 3/2 and xq = 0.3. The steric interaction, f s teric,o{d), of two monolayers through 
an oil film is obtained from Eqs. (Q) and (|^) by interchanging I and d. Thus, the total 
entropic part of the free-energy density of the lamellar phase is fthermih d) = f s teric,w(l) + 

/ steric,o\d) ■ 

To calculate the entropy of the droplet phase, we assume that the swollen micelles 
behave as hard spheres of radius R + 1/2. The free-energy density is then very well 



approximated by the Carnahan-Starling form [17 



f HS (R, p) = k B Tp (lnA 3 p - 1 + 7^—^ J (6) 

where rj = (4n/3)p(R + 1/2) 3 is the packing fraction, and A a microscopic length scale, 
which in the following is taken to be equal to the amphiphile length. Undulation modes 
on microemulsion droplets give an additional contribution to the free energy, which in the 
limit R Cq 1 has the form QTSJ 

fdrop(R,p) = p\k B T\n(R 2 ) . (7) 



Such a term becomes important for large droplet radii, and disfavors the droplet mi- 
croemulsion in this regime. The entropic part of the free-energy density of droplets then 
reads f t herm(R, p) = f HS (R,p) + f d ro P (R,p)- 

Finally, we have to calculate the contribution of the van der Waals interaction to the 
free energy. This is done within the Hamaker approach, where a pairwise form of the 
interaction is assumed || ||. The dielectric properties of the amphiphile tails (heads) 
are assumed to be identical to those of the oil (water) molecules. We subtract the free 
energy of the pure water phase, so that only the interaction between oil molecules (and 
amphiphile tails) has to be taken into account. For a molecular fluid, the interaction 
free-energy density is given by fI7 



fvdw = TP I rf3r i / rf3r 2 Ml r i - r 2lM|ri - r 2 |) (8) 
V Jv Jv 

where wir) is the interaction potential, which consists of a long-ranged attractive and a 
short-ranged repulsive part, and g(r) is the normalized density- density correlation func- 
tion, which approaches unity for large distances and vanishes for short ones. We approx- 
imate the product of these two functions by 

e 2 1 



w(r)g{r) = <f>(r) = — 



7T 2 



( r 2 + £ 2)4 r r 2 + e 2)3 



(9) 



where A is the Hamaker constant. The parameter e in Eq. (Q) determines at which length 
scale (in units of the amphiphile length) the repulsion becomes dominant. 

The calculation of the free-energy density, f intra, of the van der Waals interaction of 
the oil molecules within one aggregate is straightforward in both the droplet and lamellar 
phases (compare Ref. ||). This also holds for the inter-aggregate part of the free-energy 
density, f inter, in the lamellar phase. In the (disordered) droplet phase, the inter-aggregate 
part is taken to be 

/inter, drop = TT / <^ r W n( r ) ( 10 ) 

2 Jr>2R+l 

where Wr{t) is the interaction energy of two spheres of radius R with distance r between 
their centers. To obtain Eq. (0), we have replaced the two-particle distribution function 
for hard spheres of radius R + 1/2 by a step function |19| |. It is important to note that 
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with the choice (|9]) for the interaction potential, all interaction integrals can be carried 
out analytically. However, the results are rather lengthy and can therefore not be given 
here. 



3 Phase Diagrams 

In order to calculate phase diagrams, we add chemical potential terms to the free energy, 
and minimize with respect to I and d, or R and p, respectively. Since the functional 
form of K e ff(l;d) was determined self-consistently in Ref. [|16j, we have to use the same 
procedure here. The free energies in our model depend on the parameters k, K, Co, 
and A. The value of the Hamaker constant is well known experimentally to be A = 
(0.28T/300-K" + 0.17)10~ 20 J for the interaction of hydrocarbon across water ||. For 
T = 300K one therefore has A ~ 1/c^T. For the spontaneous curvature, we assume a 
linear temperature dependence of the form c (T) = c (Ti)(T 2 — T)/(T 2 — Ti), where T 2 
denotes the phase-inversion temperature. 

A typical sequence of phase diagrams with decreasing spontaneous curvature c is 
shown in Fig. 1 for bending rigidities of the order of lksT. The parameters are chosen 
such as to allow a comparison with the experimentally well studied oil-water-C*i2-E5 sys- 
tem. In addition to the two-phase coexistence region of lamellar and droplet phase, we 
also show the two-phase coexistence of the droplet microemulsion with a pure oil-phase 
(emulsification failure). We can follow the two-phase coexistence curves in the Li-phase 
only to a minimum water concentration of about 30%, where the oil-filled micelles become 
close-packed. The most prominent feature of the phase diagrams shown in Fig. 1 is that 
with decreasing spontaneous curvature, the lamellar phase extends further and further 
into the water apex of the Gibbs triangle. This behavior is in excellent agreement with 



experimental observations 20, 21 



Phase diagrams for larger bending rigidities are shown in Fig. 2. For n/ksT ~ 3, the 
reentrant behavior of the droplet phase disappears, and the lamellar phase is found to be 
always stable for small oil content. Finally, for n/ksT ~ 10, the lamellar phase (with low 
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oil content) can no longer be swollen to large intermembrane separations, but coexists 
with an almost pure water phase. 



4 Discussion 

The extension of the lamellar phase into the water apex of the Gibbs phase triangle 
for n/ksT ~ 1 can be understood as follows. Consider a path in the Gibbs triangle 
of Fig. lb of constant amphiphile concentration, for which the system is in the droplet 
phase for equal oil- and water concentrations. With decreasing oil content, the droplet 
size decreases, so that the bending energy cost for droplets increases, and ultimately leads 



to a transition into the lamellar phase. This is the effect considered in Ref. [Il|. With 
further decreasing oil content, the steric interaction of the monolayers across the oil films 
increases, and leads to a (reentrant) transition to the droplet phase. 

When the bending rigidity is increased to n/ksT ~ 3, the lamellar phase can be 
swollen by a very large amount by adding water. Thus, the repulsive interaction dom- 
inates in the lamellar phase for /t//c^T^3. Only when the bending rigidity is increased 
to n/ksT > 3 does the attractive interaction begin to affect the phase behavior, and un- 



binding transitions [Q can be observed (compare Refs. |52| for a discussion of unbinding 
transitions in the binary water-amphiphile system). The first-order character of the phase 
transitions for k//c^T^3 discussed above should therefore hold beyond the superposition 
approach (g). 

It is also interesting to note that the van der Waals interaction is not strong enough 
to produce a two-phase coexistence between two droplet micro emulsion, in agreement 
with recent calculations for colloidal systems This result may change when the 

deformability of the oil droplets is taken into account, which allows a larger contact area 



of two adhering droplets and thereby increases the interaction energy j23 . 

Our work can be extended in several directions; it would be interesting, for example, 
to investigate the stability of the hexagonal phase of cylindrical micelles. A more de- 
tailed comparison with experiment requires the systematic measurement of ternary phase 
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diagrams with larger bending rigidities. 
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Figure 1: Ternary composition phase diagrams for low bending rigidity and (a) T = 293K 
and (b) T = 310f^. The parameters of the model are n/ksTi = 1.0, k/k = —0.5, 
Co(Ti) = 1/6, e = 0.7, T\ = 2937^ and T 2 = 321i^. Each Gibbs triangle contains two 
coexistence regions: for higher amphiphile concentrations the lamellar phase L a coexists 
with the droplet phase L±, for lower amphiphile concentration the droplet phase coexists 
with an excess oil phase (emulsification failure). 
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Figure 2: Ternary composition phase diagrams for (a) intermediate and (b) high bending 
rigidity. The parameters are Cq 1 = 15.3, k/k = —0.5, e = 0.7, and k/IcbT = 2.8 and 
k/UbT = 9.5, respectively. 
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